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Recently the LHCb collaboration reported evidence for direct CP violation in charm decays. 
The value is sufficiently large that either substantially enhanced Standard Model contributions or 
non-Standard Model physics is required to explain it. In the latter case only a limited number 
of possibilities would be consistent with other existing flavor-changing constraints. We show that 
warped extra dimensional models that explain the quark spectrum through flavor anarchy can 
naturally give rise to contributions of the size required to explain the the LHCb result. The D 
meson asymmetry arises through a sizable CP-violating contribution to a chromomagnetic dipole 
operator. This happens naturally without introducing inconsistencies with existing constraints in 
the up quark sector. We discuss some subtleties in the loop calculation that are similar to those in 
Higgs to 77. Loop- induced dipole operators in warped scenarios and their composite analogs exhibit 
non-trivial dependence on the Higgs profile, with the contributions monotonically decreasing when 
the Higgs is pushed away from the IR brane. We show that the size of the dipole operator quickly 
saturates as the Higgs profile approaches the IR brane, implying small dependence on the precise 
details of the Higgs profile when it is quasi IR localized. We also explain why the calculation of 
the coefficient of the lowest dimension 5D operator is guaranteed to be finite. This is true not 
only in the charm sector but also with other radiative processes such as electric dipole moments, 
b — > S7, e'/ejc and /1 — > &y. We furthermore discuss the interpretation of this contribution within 
the framework of partial compositeness in four dimensions and highlight some qualitative differences 
between the generic result of composite models and that obtained for dynamics that reproduces the 
warped scenario. 



I. INTRODUCTION 



Recently the LHCb collaboration reported 3.5c evidence for a non-zero value of the difference between the time- 
integrated CP asymmetries in the decays D° — > K + K~ and D° — > tt + tt~ [1] Aacp = o-k+k- — a ir+Tr-, where 

= r( D° -> f) - T(D° -> /) 
a f ~ T ( D o _> f)+T(D° ->f) ' 1 ' 

Combined with other measurements of these CP asymmetries [2H6] , the present world average is 

Aa C p = -(0.67 ±0.16)%. (2) 

The effective weak Hamiltonian relevant for hadronic singly-Cabibbo-suppressed D decays renormalized at a scale 
m c < fj, < nib is 

^iA M d=i = % E A * E c? «? + h - c - + ■ • ■ > ( 3 ) 

q=s,d i=l,2 

where X q = V* q V uq , Q\ — (uq)v-A (qc)v-A , Q\ = {u a q/3)v-A (q~/3Ca)v-A , and ot,j3 are color indices. Dots denote 
neglected Standard Model (SM) penguin operators with tiny Wilson coefficients. In the SM, as well as within its 
minimally flavor violating extensions [7], contributions of the Hamiltonian T-lf^ =1 to Aacp are suppressed relative to 
the leading CKM terms factored out in Eq. ^ by |T46V^&|/|T4sV^ s | ~ 0.07% and are therefore expected to be small [8]. 
However, since the charm scale is not far from Aqcd> non-perturbative enhancements leading to substantially larger 
values cannot be excluded [S] (see also [TO])- 

Nonetheless, without a substantial enhancement, the SM contribution would be too small to explain current ob- 
servations. Moreover extensions of the SM generally also have difficulty accommodating the measured value without 
conflicting with existing stringent flavor-changing constraints [11] . since loop effects can produce other flavor viola- 
tions in excess of their experimental values. This however is not true for the chromomagnetic dipole operators [5J 
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due to the light quark mass suppression that is essential given the operators' helicity structure. Any new physics 
contribution to these operators can be encoded in 



\Ac\ = l 



V2 



(C 8 Q 8 + C'M + h.c. 



(4) 



where Qg = — <? s m c w^„(l + 7 5 )T a G^ I/ c and Q' a obtained from Q s with 75 — ► —7,5. The contributions of such 
operators to Aacp are given by [llj 



Aa chromo 



<-r -9Im(i?f + RI) J2 Im (C) 

i=8,8' 



where i?f denote the relevant operator hadronic matrix element ratios. Following [8] we find 

1 + rH3 



I dpi _ 167r 



1 + C 2 /3d 



12(9) for P = K(tt) , 



(5) 
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where we have used a s (m c ) ~ 0.35, C2(m c ) ~ —0.4, Ci(m c ) ~ 1.2 [13] . and = 2mp/m c (m q + m u ) ~ 3.6(1.7) for 
P = K(tt) and (7 = s(d) using inputs from |14j . We note that under the above normalization convention for the dipole 
operators in Eq. Q and using naive factorization, their matrix elements are enhanced by a factor of O(10) relative 
to the tree and penguin operators. Thus, given 0(1) Wilson coefficients of the dipole operators in Eq. Q, the scale 
required to saturate the experimental value in Eq. (2| is Ag ~ 20TeV. This is comparable to the scale A~4f ~ lOTeV 
that was found for the four fermion operators in [11|. Notice that for this calculation we have evolved the Wilson 
coefficients computed at the scale Ag down to the charm mass scale. The leading order anomalous dimension for Q 



(') 



is 7 = 28/3 in the conventions of [T5] 
where 



Consequently, the Wilson coefficients at the m c scale are (m c ) = f}C^ (A ; 



A') 1 
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II. DIPOLE OPERATOR IN RS 



One important feature of the bulk Randall- Sundrum (RS) framework [15] is that beyond explaining the electroweak 
scale, it can lead to hierarchies in quark and lepton masses without any hierarchies in fundamental flavor parame- 
ters [13 [H] • This is achieved when different flavors are localized in separate regions along the extra dimension [T7J [TH] 
through their different bulk masses. In this setup, gauge bosons, fermions, and the Higgs boson are in the bulk with 
light quarks localized in the UV while the Higgs and KK modes are localized in the IR. The assumption is that 
mixing angles arise (up to order one factors from the Yukawas) as the ratios of left-handed wave functions while 
the ratios of masses are determined by the ratios of right handed wave functions (divided by mixing angles). The 
interesting feature in the flavor sector is that although an explicit GIM mechanism is absent, flavor-changing processes 
are suppressed by wavefunctions related to quark masses and mixing angles and are therefore reasonably consistent 
with observations [19] . However, some flavor changing processes can exceed their SM values, such as Aacp which 
we consider here (for observables related to the down sector see [2"0Tf2"2"] ) . We now show that the version of RS that 
accounts for quark and lepton masses through anarchic up-type Yukawa matrices can give a sufficiently large con- 
tribution to Aacp through the chromomagnetic operators in Eq. Q. As with any model contributing through the 
chromomagnetic operator, the contribution can in principle be large enough to account for Aacp without introduc- 
ing overly large flavor-changing effects in other D meson processes [H [T^] . The important additional feature in RS 
flavor models is that all flavor-changing processes are suppressed by wavefunctions, with no additional requirement 
of near-degeneracy in masses, as there would be for example in the squark sector of a supersymmetric model. Such a 
requirement would be difficult to accommodate consistently with the requisite off-diagonal flavor entries in a model 
with only symmetries as a constraint. 

We work in a slice of AdSs spacetime whose fifth (conformal) coordinate z is bounded by two branes at zjjv ~ 
(10 19 GeV) -1 in the UV and z ~ TeV~ x in the IR. We also assume that the Higgs field, H, is a bulk field with vacuum 
expectation value (VEV) 



(8) 



FIG. 1: Leading one-loop diagram that contributes to the operator Qg from neutral Higgs exchange. Arrows denote the flow 
of charge through the diagram. Diagrams with a gluon radiated from the external quark legs are ignored as they contribute 
only to the gauge coupling renormalization. 



with v ~ 246 GeV. (3 parameterizes the localization of the Higgs VEV in the bulk. We assume for simplicity that the 
SM-like Higgs fluctuation h has the same profile as its VEV along the fifth dimension, which is a good approximation 
up to C(m^/TO^ K ) where ijjkk is the KK scale. The case j3 = corresponds to gauge-Higgs unified models (see 
e.g. [Ml). 

Since in RS the only terms that generate interactions between SU(2) doublets and singlets are the 5D Yukawa 
interactions, the Hamiltonian in Eq. Q is induced through a one-loop contribution involving exchange of KK fermions 
and the Higgs boson [19]. Some of the relevant diagrams are shown in Fig. [I] Following Refs. [2TJ [22], an explicit 
evaluation of the one-loop amplitude yields 

C 8 (m KK ) = E7» 5 2 Op , (9) 



16n 2 G F 7 



'KK 



where tokk — 2.45/ R' is the KK scale and Y5 is the 5D Yukawa coupling in appropriate units of the AdS curvature. 
Cg is obtained through the replacement Ul — > Ur and U^ R ^ denotes the 1-2 mixing angle in the left- (right-) handed 
up quark sector. Under the assumption of anarchic RS flavor parameters, the mixing angles are ~ fgi / Jqi — A c 

and U R 2 ~ fu/fc — (Tn u /m c )/X c , where A c ~ 0.23 is the sine of the Cabibbo angle, and carries an arbitrary 0(1) 
CPV phase provided by the anarchic 5D Yukawa. f x i is the value of the x l quark zero-mode profile on the IR brane. 
The contribution to the chromomagnetic operator involving the left-handed up quark, Q$, can be large in the RS 
scenario because of the sizable ratio of the first two generation wave functions, which is comparable to the relatively 
large Cabibbo angle. 

The function Or parameterizes the Higgs profile overlap with the KK state wavefunctions, and depends significantly 
on the Higgs field localization parameter (3 as we describe below. Plugging Eq. ([9| and Eq. ([6| back into Eq. (pi), 
along with RS parameters yielding the correct quark masses (see e.g. [21]) and including a running factor from 3TeV 
to the charm scale of 77 ~ 0.4, we find 

\**~\.<*<»**(%)$)'(^)'. 

which is of the right size required by Eq. ([2| . For the sake of definiteness we have assumed a maximally delocalized 
Higgs ((3 = 0), as arises in gauge-Higgs unification models. In this case the overlap function is Op ~ 0.1 for UV- 
localized first two generations. 

As we already mentioned, Op has a strong dependence on the Higgs profile. In particular as we show below for 
large values of (3 > 10 mimicking a Higgs profile very peaked towards the IR brane the overlap becomes Op ~ 0(1). 
The LHCb result is thus reproduced with a smaller 5D Yukawa of Y 5 ~ 2, similarly to generic composite models |25j . 
We study below the dependence of Op on the Higgs field localization in the bulk and we present the implications of 
A.acp within up flavor anarchic RS models as a function of the Higgs overlap function. 



A. Higgs profile dependence 

We now evaluate the RS contributions to the dipole operator for varying Higgs profiles. We describe some subtle 
features of the calculation, including the Higgs profile and cutoff dependence, and, in particular, the result with a 
properly regularized brane-localized Higgs. We address issues related to cutoff dependence, where we need to account 
for both the UV cutoff and Higgs regularization (for the brane-localized case) as discussed in [2B]- We furthermore 
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address the critical dependence on the Higgs coupling to the "wrong chirality modes" , the KK modes with opposite 
chirality to the zero-modes. We find the wrong-chirality couplings and the heavy KK modes near the cutoff play a 
critical role and furthermore one has to account for the order of removing the cutoff of the 5D effective theory and the 
Higgs width parameter. We describe here our computation method and some qualitative features of the result, which 
we develop in detail in the Appendices. We show a simple argument for the necessity of both the wrong-chirality 
Higgs coupling (which in turn helps explain the finitcncss of the result) and the critical contribution of heavy KK 
modes (those with KK number inversely proportional to the Higgs profile width) in the mass insertion perturbative 
approach. 

As in Ref. [57j , we assume the Higgs is a bulk field whose localization is controlled by the free parameter j3 defined 
in Eq. The larger /?, the more the Higgs profile is peaked towards the IR brane. We evaluate the relevant 

dipole contributions at one-loop using KK decomposition and sum the fermion KK towers up to arbitrarily high KK 
masses, thus capturing any possible sensitivity to high mass modes. Besides determining the dependence of the dipole 
contributions to the Higgs profile in the bulk, such a calculation can also serve as a regularization procedure for the 
IR localized Higgs case [27] by taking the formal limit f3 — > oo. 

For any finite f3 the KK sum is guaranteed to converge since the 5D operator g^Qa (J „G'"' Hu is finite for a bulk 
Higgs [19 . On the other hand naive 5D power counting suggest that the same operator might be logarithmically UV 
divergent when the Higgs is an IR localized field and therefore its dipole contribution, being UV sensitive, would be 
incalculable. The difference can be accounted for by the loss of momentum conservation for a brane-localized Higgs, 
which would allow two independent KK sums for fermions in the loop, thereby increasing the degree of divergence. 

A closer look at the bulk equations of motion shows that not all the fermion wave-functions are unambiguously 
defined for a delta-function Higgs. In any case, such a Higgs profile ceases to make sense for a regulated brane with 
finite thickness of order of the cutoff. Clearly the Higgs profile will not be more localized than the brane thickness. 
The IR brane Higgs can however be regularized by spreading its profile on the brane or even into the bulk [27] . A 
Higgs profile like in Eq. ^ with a large but finite f3 would provide such a regulator. Using the procedure outlined 
above and taking the large f3 limit we find that the dipole contributions converge to a finite value, from which we 
conclude that the latter are one-loop finite and thus calculable even for an IR localized Higgs. 1 We now present a 
simple argument to justify this result. 

The first important point in understanding the result is the essential insertion of the Higgs coupling to the wrong 
chirality modes. This result does not arise from a symmetry of the theory but rather follows from the particular 
structure the Yukawa coupling matrix has in the limit where the wrong chirality Yukawa is absent. 

We show that the 1PI contributions are controlled by the Yukawa interaction coupling of the heavy fermion of 
wrong chirality, opposite to that of the chiral states. Notice that a PQ symmetry does account for the need for wrong 
sign chirality in the case of neutral Higgs exchange. However, this argument does not hold for the charged Higgs, for 
which both up and down quark couplings to the Higgs field are present. Nonetheless, even when a charged Higgs is 
exchanged, the relevant diagrams vanish in absence of the wrong chirality Yukawa. 

A hint that a structure beyond PQ symmetry is responsible is seen as follows. Although 1PI diagrams with an 
up-type quark running in the loop are forbidden by the PQ symmetry for a complex Higgs, diagrams with only the 
real or imaginary Higgs component in the loop are still formally allowed. The corresponding loop amplitudes have 



the same form as in Eqs. (|A3|,(|A4j) and thus the term cx y^y}^*y}^ vanishes in the m# — > limit for both real and 
imaginary Higgs components individually. This observation suggests that in the absence of Yukawa coupling to the 
wrong chirality states both the neutral and charged Higgs 1PI diagrams vanish even without invoking a PQ symmetry, 
indicating the PQ symmetry does not suffice to understand the absence of 1PI dipole contributions arising from the 
y 11 interaction. 

We now show the particular structure of the Yukawa coupling matrices (in the absence of the wrong chirality Higgs 
coupling) leads to a cancellation of the leading Higgs dependence at each KK level. This structure is not the result 
of a symmetry of the theory but resembles the mechanism of Nelson-Barr for solving the strong CP-problem |28) . 
We focus here on the Yukawa interactions among up-type quarks and the neutral Higgs component H° = v + h, 
where v ~ 246 GeV is the Higgs VEV and h is the SM-like Higgs fluctuation. The relevant part of the Lagrangian in 



Eq. ((AT]) is 

- C D u L Y u u R h + u L M u u R + h.c. , (11) 



We note that the same regularization procedure is expected to render a finite Higgs-gluons couplings as with a bulk Higgs they arise 
from the \Y^\g^H^ HG^ U G ^ v operator, which is irrelevant in 5D. 
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where ul = (Q°, Q 1 , wl), Ur = (u°, Ql, u 1 ) and 

fy°u° o y\ 

Yu = Vi° y^ 1 I , M a = I y> m ^ I , (12) 





where we have assumed tuq = m u = m. We set now y~ = and show how the one-loop dipole contribution is 
suppressed in that case. Since m 3> v the two heavy KK states nearly maximally mix through . Alter diagonalizing 
the corresponding 2x2 block by means ol a bi-unitary transformation the heavy eigenmasses are m± = m ± y^v. 
Below we consider these two states as approximate mass eigenstates and treat their remaining mixing with the 
zero-mode perturbatively. The projections and y^° of the heavy eigenstates onto the Q° and u° zero-mode are 



,.o± _ + Vu ( , , 2A<vA i ±0 _ Uu ( i i \ , | ... 



respectively. Notice the extra relative sign between the two heavy mode projections on Q , which comes from the 
fact that one of the two unitary transformations has to involve a diagonal "phase" of 7r in order to keep the two 
eigenmasses positive. (Alternatively, one could have rearranged the states such that the mass matrix is manifestly 
positive, in which case a sign explicitly occurs in the mass eigenstate). This sign cancels against the sign of y^ 1 in 
the heavy masses m± . For each heavy eigenstate we now show that there is a cancelation at leading order in v in the 
dipole amplitude between the y^ 1 correction to the KK mass and the projections on the zero-modes. The one-loop 
dipole amplitudes is of the form [22] 

^«E^- ( 14 ) 

The leading contribution to the dipole operator contains one chirality flip and is therefore linear in the Higgs VEV 
v. One can extract this linear piece by taking one derivative of the above expression with respect to the Higgs VEV 
which yields 

j=± 

Since dm±/dv = iy^ 1 the leading order contribution to the dipole operator vanishes for each KK level. 

The above shows that the coupling of the Higgs to the wrong chirality modes is critical to the leading contribution 
to the dipole operator, and therefore the behavior of the wrong chirality modes near the IR brane plays a crucial role. 
Since the equations of motion force the wrong chirality fields to vanish (at least in the absence of a Yukawa-dependent 
delta- function source [22]), the result with a delta- function Higgs profile is ambiguous since the delta- function is 
infinite at the point where those fields vanish. This ambiguity can be resolved by the beta-function regularization 
mentioned above that gives the Higgs boson a finite thickness in the bulk. This Higgs "width" can be taken as small 
as the brane thickness which must be no greater than the UV cutoff of the theory on the IR brane. The calculation 
can then be done explicitly with five-dimensional wave-functions in the presence of the nontrivial Higgs profile. 

Alternatively, the calculation can be done with perturbative insertions proportional to the Yukawa coupling without 
solving the full 5D equations of motion. We take the latter approach here and consider the net contribution of KK 
modes up to the cutoff scale. We will see that as long as the cut off scale is much bigger than the inverse of the width, 
the result converges to a /3-indepcndcnt value, but that only heavy KK modes with masses of order the inverse Higgs 
profile width are relevant. At any large but finite /3 the Higgs overlap with fermion KK modes of high enough KK 
number starts probing the "bulky" nature of the Higgs and the KK sum converges, as dictated by 5D power counting 
for a bulk Higgs field. The finiteness of the RS contributions for any f3 appears to be consistent with the finding 
of [29]. 

The function Op in Eq. ^ collectively represents the explicit evaluation of the Higgs overlaps with the KK fermion 
wave functions as well as the summation over the fermion KK towers. The dominant diagram (shown in Fig. [IJ to 
the dipole amplitude is controlled by the Yukawa coupling to the wrong chiralities [2 2) as argued above. The overlap 
function then parametrically behaves as (see Appendix [A]) 



(16) 
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where yki corresponds to the effective Yukawa coupling between KK fcrmions of k and I fermion KK number respec- 
tively and y~ stands for the effective Yukawa coupling to the wrong chirality KK fermions. The yoo factor arises 
from the SM fermion mass replacement and we have approximated the n-th KK level masses by TO^k — nm KK- The 
numerical value of Op as a function of the Higgs localization is shown on the left-hand side plot of Fig. [2] where the 
different curves correspond to various numbers of KK modes included in the calculation. Notice that we have rescaled 
the 5D Yukawa by a factor y/1 + (3 in order to maintain the effective Yukawa finite (and (3 independent) in the large 
(3 limit [3T]. As anticipated by [37] for large enough (3 the function Op saturates and only very weakly depends on the 
precise value of (3. Practically the Higgs can be considered a brane- localized field for values as large as (3 > O(10). 

We demonstrate this result analytically in Appendix [B] for a fiat extra dimension of size L. Besides being simpler 
to analyze the flat extra dimension case provides a good description of the problem at hand here since the relevant 
dynamics corresponds to the deep UV regime of the theory and locally, near the IR brane, the details of the bulk 
geometry should have no practical impact on the result. In particular we show that for a Higgs quasi- localized near 
the IR brane with a width e <C 1/Ltokk the overlap Op receives support dominantly from a shell of ~ 1/e 2 modes 
with KK number of 0(1/ e). 

A simple argument, in the warped case, for how each KK mode contributes follows from the asymptotic form of 
the heavy KK modes near the IR brane, where they behave as sin(m n z) with z~ 1 = ke kr and m n is the KK mass of 
the n-th KK level. Let z = zexp(-e) where we have taken r — f ~ e/k ~ with z = e r /k, so that we are still 
within the region of support of the Higgs wave function. Since sin(m„z) = by the boundary condition so expanding 
the exponential and sin function we get sin(m n z) = — cos(m n z) sin(m„ze). So if m n ~ (ez) -1 we are left with the 
argument of order unity so that the heavy modes are expected to give the dominant contribution. 

Each KK mode with mass of order mKK/e contributes with a 0(e 2 ) suppression factor to the dipole operator. 
However the sum over all the relevant modes in the shell yields an e-independent result, thus explaining the saturation 
observed in Fig. [2] One sees explicitly on Fig. [5] that for a Higgs width of 1/(3 a number of 0((3) KK modes needs to 
be included in order to approach the correct answer. One can also see from the figure that the support for the KK 
fermions and Higgs overlap separates when the Higgs is peaked further in the IR. When the Higgs is pushed to the 
bulk, away from the IR brane, Op <C 1 is expected to have the following behavior: the more the Higgs is pushed onto 
the bulk, the coupling between the zero modes and the KK fcrmions becomes more suppressed since it approaches 
the limit of a flat Higgs profile for which they would be orthogonal. Furthermore, a bulk Higgs has a larger overlap 
with the elementary fermions, and the dipole operator is inversely proportional to these couplings. The coupling 
between the Higgs field and the wrong chirality KK states, that dominates the dipole contribution has non-monotonic 
but moderate dependence on /3 that is described in more details in Appendix [Bj For a maximally delocalized Higgs 
([3 = 0) and UV-localized first two generations the overlap amounts to a factor of Op ~ 0.1 while for (3 > O(10) we 
find Op ~ O(l). The region of Op for the bulk Higgs case is shown on the right plot of Fig. [21 which shows indeed 
that for low value of (3 the dipole contributions are suppressed and furthermore dominated by the first KK level. We 
find the following functional fit to the value of Op that reproduces the numerical result to better than 5% accuracy 
over the range < (3 < 50 

Op ~ F((3) = 0.71 + 0.35 tanh(0.044/3) - 0.63 exp(-0. 27(3) . (17) 

Notice that in the perturbative mass-insertion approach we have taken here, the heavy KK modes, with masses of 
the order of the inverse of the Higgs width, dominate the calculation, but the dipole operator is nonetheless controlled 
by a scale of ttjkk in the final answer. Although the answer is finite and cutoff-independent, the dipole contribution 
could nonetheless be sensitive to higher dimension operators with derivatives in the z direction, since once inserted 
into the loops we have considered those would yield an unsuppressed amplitude. However, a separation of scales 
between the inverse width and the cutoff is required to regularize the calculation consistently within the 5D effective 
theory. Thus as long as there is hierarchy between the inverse width and the cutoff, these higher order operators 
dependent on the cutoff can be neglected. 

B. Implications of Aacp and RS up flavor anarchy 

We find several interesting implications of the observed CP violation in charm decays if it arises from the RS 
conribution to the chromomagnetic dipole operator. We also outline differences between RS and generic composite 
models in 4D. First of all notice that the contribution of Q' s is negligible in anarchic scenarios due to the additional 
up-quark mass suppression which leads to a factor of (m u /m c )/A 2 ~ 1/40 suppression relative to the dominant Q§ 
contribution. Second, the central experimental value requires that the 5D Yukawa changes from ~ 6 for (3 = 
down to ~ 2 for large values of (3. The /3 = case requires 5D Yukawa value close to the perturbative bound of 
47t/v / ^Vkk — 7 for Nkk = 3 perturbative KK states. Such a large Yukawa is further motivated by the need to 
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FIG. 2: Bulk Higgs overlap function charaterizing the size of the dipole contribution at one-loop. We have rescaled the 5D 
Yukawa by a factor yjl + j3 in order to maintain the effective Yukawa finite in the large ft limit. The various lines show the 
result of Eq. (16 1 where the KK sums are performed up to 1,2,5,10,20,30,40 KK levels, from bottom to top. For values of 
j3 > O(10) (left plot) we find (1 + fi)Op ~ 0(1), the overlap is suppressed for smaller values of /3 (right plot). 



suppress the contributions to ex 21 , but is otherwise arbitrary and can potentially lead to large corrections to Higgs 
production and decay rates [55]. Such a large Yukawa would also enhance the RS contribution to direct CPV in 
neutral Kaon decays [22]. Both ek and e'/e constraints can be satisfied with tokk ^ 5 TeV and I5 ~ 5 22] . However, 
clearly such a large KK scale would induce a rather severe little hierarchy problem, so a parametric approach that 
relaxes the tension with the constraint from down type flavor violation is desired. If the bulk masses are aligned with 
the down Yukawa matrices, resulting in a model of up type quark anarchy [30] . the tension is reduced, allowing for a 
lower KK scale and a large D meson asymmetry. 

CP-violating D meson mixing is inversely proportional to the size of the 5D Yukawa in warped models |31j and 
it is interesting that the current limit on this effect is saturated for a 5D Yukawa of ~ 1.5 for the brane Higgs case 
[2~4"l |3"T] . With only mild experimental progress the LHCb collaboration would thus start probing the parameter space 
of RS with up flavor anarchy and down alignment. In other words the observation of a CPV signal in D — D mixing 
could be used to determine the degree of IR localization of the Higgs field, provided that Aacp is indeed coming from 
the mechanism discussed above. 

Note that generic composite models do not capture the distinctive feature of RS models that the composite dynamics 
is in fact chiral. That is seen by the differing shape of the KK modes for the right chirality and wrong chirality modes, 
extending the chiral structure of the zero modes. In a deconstructed version of this theory, the extra chiral nature of 
the RS dynamics means that a two-site model is not sufficient to capture the full composite dynamics. The behavior 
of the overlap function Op can be obtained in a model with at least three sites. The extra suppression resulting from 
the overlap correction Op is inherent to warped models when the Higgs is a bulk field [2TJ [22l |32] . Because the loop 
factor and the possible overlap suppressions in the RS result in Eq. ([9]) are not generic to models of four dimensional 
(4D) composite Higgs with partial compositeness, generic composite models can account for Aacp without as large 
a Yukawa |25j which makes these models in principle distinguishable from a bulk Higgs RS model, again through 
their larger contribution to CP-violating D meson mixing. However, one should keep in mind that generic partial 
composite models also predict even larger down type flavor changing neutral currents which render them somewhat 
tuned. 



III. OTHER OBSERVABLES 



Given the significant size of the chromomagnetic operator, we explore whether other large dipole contributions can 
be observable in radiative or rare semileptonic D meson decays. We will see that the RS flavor contributions although 
sizable, are suppressed compared to the long distance SM contributions that dominate for CP-preserving processes. 

The calculation can be encoded in 

^Mdipole = ^ { C 7 Q 7 + C' 7 Q' 7 ) + h.c. , (18) 

with Qj — —eQ u m c ua fJll ,(l + j5)F lJlV c, and Q 7 obtained from Q 7 with 75 — > —75. The contributions of such operators 
to radiative charm decays can be estimated in the heavy quark expansion for the charm quark. We normalize the 



radiative rate to the inclusive semileptonic rate T(c 
the known leading order result for the latter (see e.g. 
branching fraction of D° — > AT 7 as 



se + v e ) to suppress charm quark mass dependence. Using 
3]) and inputs from |14j we estimate the inclusive radiative 



B(D° -> Xj) 



r(c — > wy) 
T(c -> se+v e ) 



B(D° 



5.3 ]T l^l 2 ' 



(19) 



where at leading order T(c—> wy) = ^-G 2 F m 5 c ^2 i=7 7 , \d\ . In RS, one-loop contributions to Eq. (18) arise not only 
from diagrams similar to Fig. [T] but also from analogues where the photon is emitted from a charged Higgs running 
in the loop [3T], see Fig. [5] These yield comparable contributions and we obtain CV(tokk) = (5/2)Cs(tokk) , where 
Cs(rjiKK) is given by Eq. (J9j) . Here again dominates due to the sizable Cabibbo angle. The pairs of operators Q 7 



and Qg mix under QCD evolution, in particular C 7 ' (m c ) = r)[rjC 7 '' (itikk) + 8(77 — l)Cg ' (tukk )], where 77 is given 
in Eq. ([7]) with the replacement A 8 — > tuxk while 



a s (m b ) 


2/25 


~a a {m t )~ 


2/23 


~a s (m K K)~ 


_a s (m c )_ 




_a s (m b )_ 




a s {m t ) 



2/21 



Using the values of the above RS parameters to accommodate the Aacp measurement and ij ~ 0. 
Eq. pi) 



(20) 



we obtain from 



B{D° -> AT 7 ) 



RS 



1 x 10" 



(21) 



which is three orders of magnitude smaller than the estimates of long distance dominated SM contributions to these 
decays [31]. Nonetheless, this may be enough to probe such contributions using CP violating asymmetries [35]. In 
particular for the exclusive D — > p(u))"f final states maximal RS effects are obtained by marginalizing over the unknown 
strong phases of the interfering amplitudes, yielding 



a 



RS I 



0.03 



10- 



B{D° -> p(w) 7 ) 



1/2 



(22) 



an order of magnitude above SM expectations for these observables [3S]. Similar effects are also predicted in the 
D° —> K + K~j mode away from the <fi resonance peak. The present experimental bounds on the two dominant 
exclusive modes B(D° — > p(u)~f) < 2.4 x 10 -4 P3] are an order of magnitude above long distance estimates. With the 
large D° data samples at the LHCb and the projected Super Flavor Factories, measurements of a,s at the percent 
level should be feasible. 

The RS effects in rare semileptonic modes (like D — > ir(p)£ + £~ ) are similarly small with respect to long distance 
contributions. Compared to c -> uj, the estimate is suppressed by a factor of the fine-structure constant a. At 
leading order in inverse m c 



dT(c- 



') _ aj3 e 



ds 



24tt 



(3-#)(l-3)' 1 + - r(c->u7), 



(23) 



where 
B{D° - 
is 



(Pe 



p e +) 2 /m 2 , (3 e = yl — 4m e /s and rh e = m e /m c . Integrating over s € [4mg, 1] we obtain 
; 0.7aB(D° —> Xj). Following Eq. (21 1, the RS contribution to the above branching fraction 



B(D Q ->■ Xe + e-) 



5 x 10- 



4 M 

p°e + e 



(24) 

the long distance 
-)SM _ 10 - 6 ^ 



While the present bounds on the dominant exclusive channels are B(D° — > p°e + e ) < 10" 
contributions again dominate in the SM and are estimated to yield for example B(D° — > 

which is orders of magnitude larger than our RS short distance estimate. We further note that 7 ^ operators do not 
contribute to purely leptonic D° — > £ + £~ since the relevant hadronic matrix elements vanish by angular momentum 
conservation. 

We conclude that the best hope for identifying RS contributions to the CP violating asymmetry in nonleptonic D 
meson decays might be via CP violating asymmetries in radiative D decays |35| where RS contributions at the percent 
level are possible. Also crucial is the absence of large CP violation in D meson mixing as we have already discussed. 
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IV. CONCLUSIONS 

To summarize, we have shown that flavor-anarchic warped extra dimension models can generate sizable contributions 
to AC = 1 chromomagnetic dipole operators. We have found that with large Yukawa, varying from near-maximal with 
a maximally delocalized Higgs boson to about a Yukawa of around two for a brane-localized Higgs, such contributions 
can be large enough to induce the time-integrated CP asymmetries suggested by the recent Aacp measurements at 
LHCb and CDF. The large Yukawa coupling scenario is rather specific to the RS framework where the dipole operators 
have additional suppression over generic composite models. We also comment that since the dipole contributions 
respect SU(3) (u, d, s) flavor symmetry we expect that the individual asymmetries in D — > tttt and D — > KK should 
be of comparable size, which is a qualitative prediction of this class of models, as arc the related asymmetries for 
other excited states (for more details see [10]). 

Since the dominant contributions to CP violation in D — D mixing arise from tree level KK gluon exchange and 
are inversely proportional to the 5D Yukawa [HJ , they form diagnostics for the nature of the profile of the Higgs 
or alternatively can tell between generic composite models and the bulk Higgs version of RS as follows: the case 
where the Higgs is significantly extended into the bulk implies suppressed dipole contribution and a large Yukawa 
which consequently leads to suppressed contributions to CP violation in D — D mixing. On the other hand, in the 
brane Higgs limit or generic composite models a smaller Yukawa is required to account for Aacp and the CP violating 
contributions to D — D mixing could be large, not far from the current experimental bound. Thus, future improvement 
in the measurement of CP violation in D mixing could therefore in principle shed light on the precise nature of the 
composite theory. 

One should not ignore the fact that composite models suffer from a down quark CP problem (for recent review 
see [24] ) . An interesting possibility to consider is down alignment warped models [30] . which generically circumvent the 
CP RS Kaon problem but predict up type anarchy. This is perfectly consistent with the current experimental situation 
and the observed charm CP violation could have been predicted in such a framework well before the measurement. 
Furthermore an indirect support for down alignment and anarchy in the up sector comes from the recent Daya Bay 
and Reno results [37], which reported a relatively large value for the only (thus far) unknown neutrino mixing angle. 
A reasoning similar to that of the quark sector leads to a successful model of leptons in the context of neutrino flavor 
anarchy and furthermore requires that all mixing angles are reasonably large [38] , 
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Appendix A: One-loop dipole diagrams from vector-like quarks 

We present here a detailed calculation of the one-loop diagrams giving rise to dipole operators in theories where 
the SM chiral spectrum (Q°,u°,d°) is extended with heavy vector-like quarks (Q 1 +Q 1 l,u 1 +u 1 _,d 1 + <fl), where the 
minus subscript denotes the wrong chirality fermions. We assume only one heavy SU(2) L doublet, one up-type singlet 
and one down- type singlet. Besides its simplicity such a setup completely describes the minimal RS (composite) 
framework with one level of KK states (strong resonances). The result is then straightforwardly generalized to RS 
models where each SM fermion comes along with a tower of vector-like KK states. As we are ultimately interested in 
flavor violating dipole operators it is enough to only consider Yukawa interactions whose Lagrangian is 

-CD y° u Q Hu + y d Q°Hd + y° u 1 Q Hu 1 + y 1 u °Q 1 Hu + y^Hu 1 +h.c. 

+y° d 1 Q°Hd 1 + yfQ l Hd° + y^Q'Hd 1 + y-QlHul + y^QlHdl + h.c. 

+m Q Q 1 Q\ +m u u 1 _u 1 +m d d 1 _d 1 +h.c. (Al) 
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where H is the Higgs doublet, H = i<j2H* and mq, U: d ~ wkk 3> For the sake of definiteness we focus on one- 
loop diagrams matching to chromomatgnetic dipole operators of the form O g = g s Q Ho- )lv G fJ,L 'u where g s and G Mi , 
are the QCD gauge coupling and field strength respectively. We consider the transition u°(p) — s- H + Q°(p') + g(q), 
assuming the Higgs field carries no momentum. Lorentz symmetry and gauge invariance imply that the corresponding 
amplitude A4 can be decomposed as M. = g s T a e* fl (q)A4 11 with 

M? = A[u(p')^P R u{p)] + B[iu(p')a^q v P R u(p)] , (A2) 

where e M (<z) is the gluon polarization vector and q = p — p' ■ A and B are functions of q 2 . The Wilson coefficient is then 
C g = —B/2. There are two types of diagrams contributing to the corresponding Wilson coefficient C g , depending on 
whether the Higgs is attached to the internal or external quark line; we thus define C g — C g nt 4- C g xt . Notice that for 
both types it is sufficient to consider diagrams where the vector boson is radiated by an internal line, since diagrams 
with a vector boson attached to an external quark leg only contribute to the function A. 



1. 1PI diagrams 

We consider first 1PI diagrams where the external Higgs is attached to the internal fermion line. The relevant 
diagrams are listed in Fig|3] Notice there is no diagram with up-type quarks running in the loop thanks to a PQ-like 
symmetry preserved by the Yukawa couplings and KK mass terms, under which the Q° and the Q 1 carry charge 1, 
whereas the Higgs field and the Ql_ carry charge —1. Such a symmetry forces the up- type chiral states and heavy 
fermions of same chirality to be accompanied by the Higgs field and not its hermitian conjugate. Since there is only a 
single Higgs field, down-type Yukawa interactions break PQ and thus yield the only 1PI contributions. The amplitudes 
corresponding to the diagrams in Fig. [3] are 

iM a = / ^mvfPn] ;f + + /j 2 + _ m 4 tofr-ft + *vTP*] 

X (fc+ P 2 -^ [tgS i \k+p) 2 -m 2 Q [ ^ PRHp) k*—^> (A3) 
f d 4 k . „i iilk+f' + md) i(ti + p' + m d ) 

iMb = / ^*VW** p dwT% -i T^fc S 

J (27r) 4 (k+p') z — m d (k+p') z — m d 

x[iu?*P L + iy?PR] 1* + ? 2 + mQ ] [^PbHp)^^- , (A4) 
[k + p) z — ttiq k z — m H 

where Pl(r) is the left-handed (right-handed) chirality projector. We kept the loop momentum integrals in four 
dimensions as we are only interested in finite dipole contribution. After integrating over the loop momentum we find 
327r 2 C„ int = y^yf^yfil't + 4 int ) + vfvTyfiJ^ + ^ nt ), where 



dxdydz 



l-3(x + y) zm 2 Q + q 2 z(x + y)(l-x-y) 



A A 2 



Jj nt = m Q m d J dxdydz {x+ ^J rZ) , (A5) 

with A = mQ(y+z)+xm 2 i +rnj 1 (l — x — y—z)—q 2 z(x+y). If, and J\, are obtained from I a and J a with the replacements 
{Xyirid} -H- {z,toq}. The remaining Feynman integrals are over < x < 1, < y < I — x and < z < 1 — x — y. 
Integrating over the Feynman parameters and taking the limit of mjj, q 2 <C ttlq d yield 



r int ^ t int ^ sr\ H 1 . j int ^ _ "o j lnt ^ ■ ' — o ( AR\ 

a - » ~ \m%. d ) ' " 2m 2 (r 2 -l) 2 ' Jb ~ 2m 2 (r 2 - l) 2 ' [ ° j 

with r = mQ/m d . We see that 1™^ vanish in the limit where m# — > (see also [39 ). In that limit we thus find 

C int = - y™ ly<i y d +o( m ^ ) (A7) 
9 647r 2 m Q m u \ m Q,d) 
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FIG. 3: One-loop 1PI diagrams matching to chromomagnetic dipole operators. 



2. Non-lPI diagrams 

There are additional non-lPI diagrams contributing to the operator O g below the heavy scale tokk- They correspond 
to diagrams with a Higgs vertex on the external quark lines and an intermediate heavy quark propagator, as shown 
in Fig|4j The corresponding amplitudes are 

x [iy%PL + iy^Pn] % + "J [ivi°PnHp) , (A8) 



J (27t) 4 p z — m z (k+p') z — m 



Q 



x [ig s T°f] + ^ mQ ] [iyl°P R H P ) r^-^ . (A9) 
(k + p) z — rriQ k z — m H 

Notice that contributions cx y^yudVud vanish when the external quark leg is on-shell because the corresponding 
diagrams are proportional to the external quark equation of motion pu(p) or u{p')p' . Hence the remaining non-lPI 
contribution to the Wilson coefficient is 327r 2 C g cxt = yfyu*yf{J^e + J uj) + V^Vd*vf J d£' where 

j-=^/W^, j~;=^/W^, (mo) 

m Q J A u ' J m u J A Q 

and J^f = Ju X e(u — > d) with Ax = m x (x + y) + m 2 H {l — x — y) — xyq 2 , X = Q,u, d. Integrating over the Feynman 
parameters and taking the limit m 2 , q 2 <C toq u d , we find 

T CX t r^j T CX ^ ^ \ J cx t ^ \ ( A 1 1 \ 

' J 2m Q m u 2m Q m d 
The non 1PI diagrams thus yield the following contribution to the chromomagnetic dipole operators 

^ext _ x ilu Uu ilu i x Uu y d y d /i, \ 

a S2ir z mQm u d4tt z mQm d 

3. Matching to RS models 

We finally comment on RS models, which contain several towers of vector-like KK states: Q n , u 11 and d n of masses 
tUq^ u d ~ mriKK. The above result is easily generalizable to those models. Upon resuming the doublet and singlet 

KK towers the Wilson coefficient is C g = ^ n m C g n ' m ^ where 

^j(n,m) &u y d y d _j_ yu yu yu (A13) 

32ir 2 mQ^m d m ^ 32ir 2 m^rriu^ 

Since the two terms above involve a priori independent phases they cannot be combined. We chose to focus on the 
y^ contribution in the main text since flavor anarchy is only assumed in the up-sector. 
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FIG. 4: One-loop non-lPI diagrams matching to chromomagnetic dipole operators. 




FIG. 5: Additional one-loop 1PI diagram inducing magnetic dipole operators. 



4. Electromagnetic dipole contributions 



When considering QED magnetic dipole operators like 1 = eQ°H(T^F^u there is another 1PI diagram where 
the photon is radiated by the (charged) Higgs line, as in Fig. [5j The corresponding amplitude is 



d4k u{p')\ivfP R ^ + rnd) ^- lu 



2tt) 4 

x[iyl°P R ]u(p)[ieQ H 



k 2 — m 2 . 



m*PL+iv?p R ] 



i 2 e* ■ (2k-p-p') 



[{k-p) 2 -m 2 H ][(k-p>) 2 -m 2 H ) 



(A14) 



which yields a contribution to the 1 Wilson coefficient of 327r 2 C™ t D y^ 1 2/ d 1 * ^ c lnt + V^Vd *Vd°Jc^ t where 



- int 



dxdydz 



3(y + z)-2 q 2 yz(\-y~z) 



A- 



A 2 



, J c mt = vriQ-md J dxdydz 



(y + z-l) 



A 2 



(A15) 



with A 7 = m d {\ — x — y — z) + xttCq + m 2 H {y + z) — q 2 yz. After Feynman parameter integration we find in the limit 



m 2 H ,q 2 < m 2 Q d 



O 



Q,d, 



2niQm d 



(A16) 



Here again, I c lnt vanishes in the m# — > limit. Combining the amplitude from the diagram in Fig [5] with those of the 
diagrams in Figs. [3]|4j where the gluon is replaced by a photon, yields 



1 Vuv^yf , 1 VuVu*Vu 



oi„,-*„,ro 



1927T 2 niQirid 48ir 2 rriQm u 



(A17) 



13 



Appendix B: A toy flat extra dimension model 



We prove here why the overlap function Op defined in Eq. ( 16 ) is finite and does not depend on the Higgs width when 
the latter is quasi-localized on the IR. We show this result in the simple case of a flat extra dimension. Nonetheless 
we stress again that the flat case should provide a reasonable approximation to the warped case as far as the width 
of the Higgs profile is small so that deviations from the flat geometry are subdominant. Let us consider a flat extra 
dimension of size L where the wrong chirality modes are given by a set of sine functions, whereas the right modes 
correspond to towers of cosine functions. 

The Higgs profile is taken to be peaked towards the IR like e~ z ^ eL , where z is the extra dimension coordinate. We 
focus on the contribution to the dipole operator in the limit where e< 1. Since the cosine wave-functions of the right 
chirality modes are approximately flat in the region where the Higgs profile is peaking, their couplings to the zero 
mode have very weak dependence on the Higgs profile in the small e limit. The Higgs overlap function Op is then 
proportional to 

e^Ol"* oc e- 1 V / dz sin (nirz / L) sm(mirz/ L) e- z / tL , (Bl) 

p mn In 

n,m=l JIJ 

where we divide both sides by e in order to maintain a finite value for the Yukawa overlap integral, equivalent to 
rescaling the 5D Yukawa by y/1 + /3 for the warped case with the polynomial Higgs profile in Eq. Q. For a high 
enough KK level, then, one can replace the sum over n and m by a continuous integral. Moving to radial coordinates, 
defining r 2 = n 2 + m 2 and tan# = n/m, the resulting integral can be solved analytically. To understand why the sum 
is finite and independent of e let us first integrate over z and leaving an r-dependent integral of (to leading order 
in e): 

c-Woc r rdr (OT } 2 (B2) 

P JV2 (l + e 2 ir 2 r 2 )Vl + 2e 2 TT 2 r 2 V ; 

The change of variable x = enr explicitly shows that the result is indeed e-independent as anticipated. Furthermore 



one can see from Eq. (B2| that the dominant contribution to the KK sums comes from heavy modes with r ~ 1/e, 
however the width of the distribution is such that 1/e 2 modes contribute to the integral, yielding an e independent 



result. Eq. (B2| also shows that for any small but finite e the ultra heavy modes with r ^ 1/e give a negligible 
contribution to the integral, because such modes rapidly oscillate within the Higgs profile's width, leading to a single 
KK sum and thus a highly suppressed contribution. This fact also explains why the total dipole contribution is 
guaranteed to be finite. 
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